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^li 1 Introduction 



Abstract. We discuss the calculation of the leading hadronic vacuum polarization in lattice QCD. Ex- 
ploiting the excellent quality of the compiled experimental data for the e^e~ — >■ hadrons cross-section, we 
predict the outcome of large-volume lattice calculations at the physical pion mass, and design computa- 
tional strategies for the lattice to have an impact on important phenomenological quantities such as the 
leading hadronic contribution to {g ~ 2)^ and the running of the electromagnetic coupling constant. First, 
the R(s) ratio can be calculated directly on the lattice in the threshold region, and we provide the formulae 
to do so with twisted boundary conditions. Second, the current correlator projected onto zero spatial mo- 
mentum, in a Euclidean time interval where it can be calculated accurately, provides a potentially critical 
test of the experimental R{s) ratio in the region that is most relevant for {g — 2)^. This observation can 
also be turned around: the vector correlator at intermediate distances can be used to determine the lattice 
spacing in fm, and we make a concrete proposal in this direction. Finally, we quantify the finite-size effects 
on the current correlator coming from low-energy two-pion states and provide a general parametrization 
of the vacuum polarization on the torus. 
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In quantum field theory, the information encoded in the 
correlation functions of conserved currents has important 
phenomenological applications. The correlation function 
of the electromagnetic current, in particular, quantifies the 
polarization of the 'vacuum' by virtual particles induced 
by the passage of a photon. This virtuality-dependent vac- 
uum polarization n{Q^) affects the propagation of the 
photon and has physically observable consequences, see 
for instance [1 . One of them is the running of the fine 
structure constant a{Q'^), which now depe nds on the four- 
momentum squared of the photon via Eq. ( |11[ ) below. An- 
other consequence of the polarization of the vacuum is a 
contribution of all virtual particles to the magnetic mo- 
ments of leptons, which can be measured to very high pre- 
cision in the case of the electron and the muon (see [2] for 
a review of the subject). We will discuss this application 
extensively below. 

One important contribution to the vacuum polariza- 
tion comes from QCD. At high virtuality Q^, this con- 
tribution to n{Q^) is calculable in perturbation theory, 
due to the asymptotic freedom property of QCD. Below 
a scale of a few GeV, the vacuum polarization receives 
large non-perturbative contributions, making it inacces- 
sible to known analytic methods. From here on, we will 
focus exclusively on the QCD contribution to n{Q^), but 
will keep using the same symbol. In particle phenomenol- 
ogy, it has been customary to extract the low-Q^ part 



of the function iT(Q^) from experiments via a dispersion 
relation (Eq. [16] below) . However it is also possible to di- 
rectly calculate 77 (Q^) from first principles using numer- 
ical lattice QCD methods (SIIHISIIH] , roughly for a range 



of momenta O.lGeV^ < Q < 4.0GeV^. What hmits the 
upper end of the momentum range is the size of the lat- 
tice spacing a, whose inverse provides a momentum cutoff 
^ ir/a. On the lower end, it is the discreteness of the 
available momenta in a finite volume (IQminl = -^ on a 
torus of dimension L) that limits the accessible Q^ values. 
We note that current correlators involving heavy flavors 
of quarks have been used on the lattice for other purposes, 
namely determining the charm quark mass and the QCD 
coupling constant [7]. 

To obtain the leading hadronic contribution to the 
anomalous magnetic moment of the muon a^^*^, the imag- 
inary part of the vacuum polarization at timelike mo- 
menta is folded with an analytically known QED ker- 
nel [5] which involves only the scale m^. The same quan- 
tity a^ can also be expressed as an integral over space- 
like momenta ^9]. An important observation of Blum [3] 
was that this opens the possibility to evaluate it in the Eu- 
clideanized theory, which can be simulated non-perturbatively 
by Monte-Carlo methods [3,4,5,6 . Because the muon mass 
/i w 105MeV is small on hadronic scales, the contribution 
to a^ is dominated by the region of small Q" or, alter- 
natively, by long distance contributions. How exactly this 
affects the calculation is one topic we will address. 
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Since the vector spectral function is already extremely 
well known from particle physics experiments, only highly 
accurate lattice predictions will have a useful impact on 
phenomenology. By the same token, the knowledge of the 
vector spectral function allows us to design useful observ- 
ables which are both accurately computable on the lattice 
and critically challenge the experimental measurements 
that are relevant to the determination of Aa{MV) and 

HLO 

"a* • 

One of our goals will be to understand some of the 

systematic errors that may affect the direct lattice cal- 
culation of the vacuum polarization. A second goal is to 
propose other promising computational strategies, in in- 
creasing level of 'ambition', which we believe could have 
a phenomenological impact. These strategies represent a 
compromise between computational feasability and phe- 
nomenological relevance. The third topic is our proposal 
of a new reference scale tq which is both accurately cal- 
culable on the lattice and extractable from experimental 
data with negligible model dependence. The time-scale tq 
is defined from the isospin vector current correlator. 

In section (pi) , after a review of the basic relations be- 
tween the relevant observables, we provide a comparison of 
the vacuum polarization calculated on the lattice and the 
vacuum polarization obtained from the phenomenologi- 
cal R ratio by employing the dispersion relation ( 10 1 . We 



observe that even state-of-the-art lattice calculations still 
yield a vacuum polarization about a factor of two smaller 
than phenomenology indicates in a wide range of mo- 
menta. In section ([3]) , we provide the general parametriza- 
tion of the polarization tensor 11 ^^(Q) compatible with 
the symmetries of the torus and the conservation of the 
electromagnetic current; this allows one to perform sev- 
eral tests of finite-volume effects. In section Q we will 
switch to the mixed representation {t,q), where t is Eu- 
clidean time, which facilitates the interpretation of the 
vector correlator in terms of physical states. At large times 
t, the current correlator is exponentially dominated by the 
low-lying states. These are two-pion states on the torus 
and due to the sparsity of momenta available to them, 
their contribution is affected by 0(1) finite volume effects. 
Based on this understanding, we will consider three ways 
in which the lattice could have an impact on the determi- 
nation of Aa{Q'^ = M|) and aJJ^°. Section ^ describes 
the new reference scale tq. In section (|6|, we present a 
technical generalization of the Liischer formula in the vec- 
tor channel for twisted boundary conditions (see also [TUl 
lllj in the scalar channel), which are expected to help in 
the determination of the timelike pion form factor. The 
concluding section contains a summary of our findings and 
proposed computational strategies. 



2 Preliminaries and status of vacuum 
polarization determinations 

We begin this section by introducing the relevant quan- 
tities and reviewing the most important relations among 
them. This will allow us to perform a comparison between 



lattice experimental data. Our Minkowski metric conven- 
tion is (H ). Euclidean momenta are denoted by a 

capital letter, Minkowski momenta by a small letter. 
The electromagnetic current 



JfT = ^uif^u - ldjf,d - is7^s + . . 



(1) 



is the central operator of interest in this paper. The cor- 
responding spectral function is defined as 



1 

2^ 



p^,{k) EE - / d4xe*'=-(0|[j7-(x),jr(0)]|0). (2) 



Due to current conservation and Lorentz invariance, the 
tensor structure of p^^ is 



Pp.u{k) = {k^k^ - g^^yk^) ■ pik^)- 



(3) 



The spectral density p is non-negative. In the free theory 
for massless quarks of charges Q/, it is given by a step 
function, 

p{s) — 9{s) (free massless quarks). (4) 

More generally, p{s) is related to experimental observables 
by the optical theorem, 

, , R{s) ^, , o-(e+e~ — >■ hadrons) ,_, 



127r2 



47ra(s)2/(3s) 



The denominator is the treelevel cross-section e+e~ — )■ 
p^ H~ in the limit s ^ m^, and we have neglected QED 
corrections. At low energies, the spectral density is given 
by the pion form factor [2], 



This relation holds near threshold, 27Tijr < ^/s < 3r?z^, and 
even up to Amj^ if the electromagnetic current is replaced 
by the isospin current in the definition of p{s). In |12j . a 
formula relating the pion form factor to a finite-volume 
matrix element calculable in lattice QCD was derived, 

iF.(i.)P = (#'(,) +fc^)^iA,r (7) 

Here E equals the invariant mass of the two pions, k is 
related to E via E = 2^vn\ + fc^, ^i is the scattering 
phase shift in the p-wave, isospin / = 1 channel and A^ is 
a vector-current matrix element between the vacuum and 
a unit-norm two-pion state j?/'") of energy E on the torus. 
Finally, g = |^ and (/> is a known kinematic function |13) . 
The scattering phase (5i(fc) can be extracted (see [I41IT5] 
and Refs. therein) from the finite- volume spectrum using 
the Liischer formula [THlll3| . 

In Euclidean space, the natural object is the polariza- 
tion tensor 



i7^.(Q)^ d'x^Q--{j^{x)j,m, 



(8) 



David Bernecker, Harvey B. Meyer: Vector Correlators in Lattice QCD: methods and applications 



and 0(4) invariance and current conservation imply the 
tensor structure 



Ht^AQ) = {Qi^Qu - s^,Q^)n{Q^). 



(9) 



The function II{Q^) can be calculated in lattice QCD [3, 
|4l[5l[6]. The leading hadronic contribution to the vacuum 
polarization e^77(Q^) in the spacelike domain can be ex- 
pressed through the spectral function via a once-subtracted 
dispersion relation, 



/■OO 

iT(Q2) _ 77(0) = q2 / j^ 
Jo 



P(^) 



S(S + Q2)- 



(10) 



An important physical application of the vacuum polar- 
ization is the running of the electromagnetic coupling, 



AaiQ"^) = 47raRe [n{Q^) - 77(0 



a(0), (11) 



In particular the value of the coupling at the scale Q^ — 
A7^ is a precision observable that, combined with the 
Fermi constant, the Z boson mass, the quark masses and 
the Higgs mass, lead to a prediction for the Weinberg an- 
gle, which confronted with its direct measurement leads 
to an upper bound on the Standard Model Higgs boson 
mass (see for instance [T]). 

The 0(0^) hadronic contribution to the muon anoma- 
lous magnetic moment can be expressed in terms oi II{Q^) 
as 



HLO 



(-) J dQ^KE{Q^)n{Q'), (12) 



with the kernel given bjr] [3] 



Ke{s) 



Z{s) = - 



-^■s-Z(sf- 

s - ^/P + 4s 
2s 



1 - sZ{s) 
l + sZ{sy ' 

s 



(13) 

(14) 
(15) 



Eq. (12 1 is used to obtain a^^'~' based on the vacuum 



polarization computed on the lattice. 



2.1 Confronting the vacuum polarization from the 
lattice and the R(s) ratio via the dispersion relation 

The vacuum polarization can be easily calculated from the 
7?(s) ratio by using the optical theorem (|5| and the disper- 
sion relation (10 1. In the Euclidean domain, the relation 
reads 



n{Q') = 



Q' 



R{s) 



"^'5(5 + Q2)- 



(16) 



^ Our function Ke matches the function / introduced in [3] . 





a 


M,/GeV 


n/GeY 





655.5 


0.7819 


0.0358 


1 


8.5 


0.7650 


0.130 


2 


11.5 


0.7820 


0.00829 


3 


50.0 


1.0195 


0.00426 



Table 1: Parameters used in the parametrization (18 1 of 
the R{s) ratio. 




1.0 1.5 2.0 

s [ GeV"] 

Fig. 1: Parametrization of the R{s) ratio. 



By taking the phenomenological determination of R{s), 
the integral can be evaluated and a comparison to the vac- 
uum polarization calculated on the lattice is possible. For 
this purpose we parametrize the R{s) ratio using Breit- 
Wigner curves of the form 



fiV^) 



cr^ 



4(^-M)2 + r2 



(17) 



where the parameters C, M, 7^ are used to match the height 
of the resonance to the experimental data compiled by the 
Particle Data Group (PDG) ,17 . Altogether our parametriza- 
tion reads, for s in units of GeV^, 



7?(s) = e{y/s - 2m^±) 6'(4.4m^± - yG) 



(18) 



3/2 



(0.6473 + /o(y5)) 
+f3iVs) + 3((|)2 + (i)2 + (1)2)0(75 - M3) 



To the fi correspond the parameters {Ci, Mi, Fi} listed in 
Tab. (0. 

A comparison between the parametrization and the ex- 
perimental data compiled by the PDG can be seen in fig- 
ure [Tl The vacuum polarization resulting from integrating 
our parametrization of R{s) in Eq. ( 16 1 is shown in figure 



(121). With the calculated vacuum polarization it is also pos- 
sible to calculate the hadronic contribution to the anoma- 
lous magnetic moment of the muon a^^'~' by using formula 



(12 1. Using our simple parametrization of the R{s) ratio. 

It is close enough 
(690.3 ±5.3) -10-1° [19] 



we obtaine the value a^^° = 672- 10"^° 
to the reference value of a^^^ 
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its anomalous magnetic moment is given by 



,HLO 



G [GeV 



Fig. 2: Comparison of the vacuum polarization calculated 
using our parametrization (181 of the R(^s) ratio and lat- 
tice data from [THUS]- 



(see the review [2] and the recent pUll^ for the latest 
evaluations) for the purpose of this article to justify the 
simple functional form we have used. 

The phenomenological hadronic vacuum polarization 
is compared to lattice QCD data generated by Delia Morte 
et al. [TSlIeFJin Fig. ^. The data is generated from A^f = 2 
gauge fielcfconfigurations, i.e. they only include the sea 
quark effects of the u, d quarks. In the electromagnetic cur- 
rent, the contributions of the up, down and strange quarks 
were taken into account. Results at different light-quark 
masses are displayed in Fig. (pi), corresponding to 'pion' 
masses m^r down to SOOMeV. The E lattices are 32'^ x 64 
and the F lattices are 48'^ x 96, so that the smallest value of 
rn-i^L is 4.7. Only the Wick-connected contributions were 
included in the calculation. Twisted boundary conditions 
were used to obtain the vacuum polarization at a denser 
set of momenta. For more details, we refer the reader to 
the original publications [T51I5]. 

The subtracted vacuum polarization n{Q^) calculated 
on the lattice lies about a factor two below the curve ob- 
tained using the experimental R{s) ratio. As the quark 
mass is lowered, the curves (very slowly) approach the 
phenomenological curve. So in spite of the excellent sta- 
tistical quality of the lattice data, the large volumes used, 
and the dense set of momenta, a large difference remains 
between the lattice and the phenomenological curve. And 
this occurs in a region (say, 0.4GeV^ < Q^ < 2.0GeV^) 
where we would expect cutoff effects to be small. We note 
however that in view of the form of kernel K{s), this re- 
gion makes only a modest contribution to a^^'~' , and that 
the smallness of a^^'~' calculated around m„ = SOOMeV is 
instead related to the behavior of IliCp) very close to the 
origin (barely visible in Fig. pi). In fact, in the limit where 
the lepton mass goes to zero, the hadronic contribution to 



^ A more extensive set of these lattice results was reviewed by 
J. Zanotti in a plenary session of the Hadron 2011 conference, 
Munich, 13-17 June 2011. 



hm ^ 

m,,->-0 raj 



n{s) 



s=t) 



(19) 



Feng et al. [5^ followed an approach based on the idea 
that using the variable Q"^ /rrii for the horizontal axis, 
where nip is the quark-mass dependent p-meson mass, 
leads to an approximate scaling at small Q^, in the sense 
that the curves corresponding to different quark masses 
would approximately lie on top of eachother. This idea 
may prove helpful in carrying out the chiral extrapola- 
tion, but it only partly explains the difference between 
the curves in Fig. (pi). 

The different channels contributing to R{s) and their 
relative importance for the evaluation of the vacuum po- 
larization and a^^'^ have been described in the literature, 
^ M 

Refs. [22, 2D; 21, 2J among others. For the reader's conve- 
nience we summarize some of the facts known from e^e~ 
annihilation experiments that are relevant to our discus- 
sion. The total two-pion contribution represents roughly 
75% of ajj""*^, and it is itself dominated by the p(770) me- 
son contribution: the part of it that comes from below 
O.SGeV amounts to 8.4% [23], while the two-pion contri- 
bution to a^^'^ from 0.1 < s < 0.85GeV2 is 69% [H 
I25| . See in particular the accurate initial-state radiation 
(ISR) measurements in this region [2^1^ : the world's 
data is nicely summarized in Fig. (3) of Ref. [21]. The 
point has also been made [27] that if the energy region 
below SOOMeV was alone responsible for the current dis- 
crepancy between the Standard Model prediction and the 
direct measurement of a^^^ ^ it would require a 52% in- 
crease of the e+e" -^ tt+tt^ cross section, and it would 
then also lead to a lowered Higgs mass upper-bound of 
143GeV at 95% confidence level. Such an increase seems 
very unlikely in view of the quoted experimental errors. 

For comparison, the entire three-pion contribution to 
a^LO from threshold to 1.8GeV is about 6.8% [2ll|20] 
(see Fig. (13) of [H]). From threshold to 660MeV, Hagi- 
wara et al. [2T] estimate it in Chiral Perturbation Theory 
to be negligible, 0.001%. The four-pion contribution has 
also been determined. From 0.305 to 1.8GeV, the 27r+27r~ 
channel yields 2.0%, while 7r+7r-27rO yields 2.7% [2l]. 

We note that a lattice calculation containing just the 
Wick-connected diagrams of the up and down quarks amounts 
to working with the isospin current. Such a calculation 
can perfectly well be compared to the experimental data 
if one selects unit isospin final hadronic states. In particu- 
lar, the w and cj) resonances should not be included in the 
dispersion integral in such a comparison. See also [25] for 
an analysis of the Wick-disconnected diagrams in chiral 
perturbation theory. 



3 The vacuum polarization on the torus 

As in all numerical lattice QCD calculations, the vacuum 
polarization is evaluated in finite volume. We will assume 
that the boundary conditions are periodic in all directions. 
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except that the fermions have antiperiodic boundary con- 
ditions in the time direction (the system is thus at a fi- 
nite, albeit low temperature). The question then arises, 
how large the finite-size effects are, and what their para- 
metric dependence is. A first step towards answering this 
question is to provide a general parametrization of the 
polarization tensor on the torus of dimensions f3 x L^. 

The correlation function of a (not necessarily conserved) 
current forms a rank-two symmetric tensor which can be 
uniquely decomposed into a traceless and a scalar compo- 
nent, 



dimensions 1, 


1, 2, 3 and 3). We then have the decompo- 


sit ions 






T^''^=Ti®T2, (27) 




rp^ ^Ai®E. (28) 


Thus 




Ht^Aq) = 


- n^^ (?) + n^Aq) + 7Tj;(g) + i7j2(g), (29) 


^^.Aq) = 


-S^^.tl^Aq)- (30) 



7T^,(g) = 7T^,((Z)+i7^,(g), 



(20) 



fi,M) = 5^un{q), El=on,,{<i) - o. (21) 



In infinite space 
implies that 



the 0(4) symmetry of the theory 



n^Aq) = {q^qu - ^s^,q^)f{q^), (22) 



n^uiq) = &t,^q^g{q^). 



(23) 



Imposing the conservation equation q^n^i,{q) = 0, we ob- 
tain the condition 



9{q') = ~f{q')- 



Returning to Eq. pO] ), one then arrives at ([9]), 
n^u{q) = [qt^qv ~ 5^,^q^)f{q^), 



(24) 



(25) 



and f{q^) can be identified with the vacuum polarization, 
usually notated n{q^). 

On a four-dimensional torus of dimensions i^, the rel- 
evant symmetry group is the hypercubic group 11(4). We 
will follow the notation of Ref. [S^. The 20 inequiva- 
lent irreducible representations of H(4) are denoted by 
Tj: , where I is the dimension of the representation and 
k = 1,2,... distinguishes inequivalent representations of 
the same dimension. There are four one-dimensional repre- 
sentations, two of dimension two; four of dimension three, 
four, and six; and two of dimension eight. The defining 
representation is labeled as rj; . The currents j^{x) be- 
long to this representation. Their direct product can be 
decomposed according to [29] 



.(4) 



.(4) 



.(1) 



.(3) 



.(6) ^ .^(6) 



(26) 



Because the vector correlator is symmetric in the space- 
time indices fi and v, the antisymmetric representation 
r} will play no role in the following. 

Very often the space on which QCD is simulated is 
a four-dimensional torus of dimensions /3 x L^. Then the 
symmetry group H(4) is further reduced to Z{2) x H(3), 
where Z{2) corresponds to Euclidean-time reversal and 
H(3) is the symmetry group of the cube. Two irreducible 
representations of H(4) appearing in Eq. ( 26 ) further break 
up into smaller irreducible representations of H(3). The 
latter are denoted Ai, A2, E, Ti and T2 (respectively of 



In matrix notation, we now have (all matrices are sym- 
metric) 



'n, 



00 







^uJl?) 



n^Aq) = 



nliiq) = 



3 -'^00 



At 



^ijAi 



(31) 



'0 
TTfi 



TTf, 

TjE 

-'-'33' 



'22 



'0 /7^i 77 -^1 77-^1 
^ -'-'01 -'^02 -'^03 



'01 





/O 



nliiq) 



nl^ 







\ 



V 



7T_ 

fjT2 

^^23 




\ 



3 ^-'00 / 

EL^^ =0,(32) 



(33) 



(34) 



Without using the conservation of the current, 7T^t^((7) is 
thus characterized by five independent functions belonging 
to various irreducible representations of the group 11(3). 

The Ward identities still read q^in^ij{q) = on the 
torus. Writing them out ioi v — {) and 1, 



9o 



9i 



'n 



Ai 






n 



Ai\ _ 



V^3 fv' 



Oj 



\n^^') = -q^n'^l-q2nl 



q^nli 



(35) 
.(36) 



There are two further equations obtained from Eq. ( |36[ ) by 
cyclic permutation of the indices (1,2,3). The three repre- 
sentations parametrizing the diagonal components can be 
related to the two off-diagonal representations Ti and T2. 
For a generic momentum with non-vanishing components 
in each direction, we have the relations 



47Toi,^ {q) - ^ ^^^ * ^ ^^1 '■> '^^ 



i:mi-i:ni}i?f-3^ (37) 






qj 



E^. 



Ti f qj 

Oj 



90^ 

qj' 
90 qj ' 



(38) 



If one or more components of the momentum vanish, one 
has to inspect the Ward identities on a case-by-case basis. 
The component in the representation E is obtained by tak- 
ing the linear combination (72 -Eq. (36) — (7iEq. (36lqj_i.g2^q3 
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The result is 



series and find for the other coefhcients 



37Tfi(g) 



91 °^ 



92 



90 ^Ti 
93 



(39) 



-77 



12 



f ^ _ 2^) + 7Tg f ^ - 2^) + 77,^^ f *^ + ^V 
^92 91-^ ^93 91^ ■ ^92 93^ 



We note that for many actions at finite lattice spacing, the 
Ward identities read 9^77^^ = with q^ = 2sin{aq^/2). 
In equations ( 35 ) and ( 36 1 , and in the equations in the 
next subsection, q^ should then be replaced by q^. 



3.1 Parametrization of the polarization tensor at small 
momenta 



Of particular relevance is the behavior of the polarization 
tensor at very small momenta, since in infinite volume 
77(0) must be subtracted from n{Q^) in order to remove 
a logarithmic divergence. On the torus we Taylor-expand 
the different contributions to the polarization tensor up 
to quartic order, including all the polynomials compatible 
with the respective cubic representations. 



3„2\ 



n^l (9) = 9091 [M + Bt^ (9? - y 



^Cti9o 



^S(9) 



= 9192 [At2 + Br^ql + Cr^q^ 

BEqoiQi^W 



VDtM 
1 „2^ 



-'^11(9) — ^-E(,9i ~ 3?' ) I ^-e-yov^l 3"^ 

+CE{qt-\Y.,1^~W(^l-y')) 
+DEq\ql-\q'), 

77it (9) = Ml + Bq^ + C'9o' + ^90?' 

h7?(g^)^ + F(E,g|-|(g^) 



,2^2^ 



n'''{q)=Aql + Bq^ + Cq 



- Dqlq^ 



(40) 
(41) 
(42) 

(43) 
(44) 



-E{ef^HT,<i^-U<i 



r,2\2\ 



It should be noted that since the momentum variable 9 
assumes only discrete values on the torus, one can always 
represent a function of 9 as a polynomial. Also, in infinite 
volume, the spectral representation ( 10 ) shows that for a 



spectral function admitting a mass gap, all the derivatives 
of the vacuum polarization exist and are finite at the ori- 
gin, and that the vacuum polarization can be represented 
as a polynomial in Q^ locally around the origin. There- 
fore one expects those coefficients in Eq. ( 40 - 44 ) that are 



commensurate with the infinite-volume tensor structure 
(|9| to smoothly tend to a Taylor coefficient of n{Q^) in 
the infinite- volume limit, whereas the others must vanish 
in the same limit. Before considering this limit however, 
we work out the consequences of the conservation of the 
current on the coefficients. 



From the Ward identity (35), we can express all five 



series expansions in terms of two of them, by relating the 
coefficients. We choose 77^^ and 77-^^ to be independent 



Ae 


= At„ 


(45) 


Be 


— Bt2 — Bti , 


(46) 


Ce 


= 2Dt,, 


(47) 


De 


= Ct2 - 7-DT2; 


(48) 


A^ 


~\At,, 


(49) 


B = 


KAt, - IAtJ, 


(50) 


C^ 


-\C!ti, 


(51) 


D = 


ii^Ti + 15-Bti - Dti - ^Bt^), 


(52) 


E ^ 


i{Dti + T05DT2 ~ ^Ct2), 


(53) 


F = 


j{Bti - 3DT2); 


(54) 


A^~ 


~|^T„ 


(55) 


B = - 


-1{At,+2At2), 


(56) 


C = - 




(57) 


b = - 


-\{-IBt,+3Dt,+Ct,+2Bt2), 


(58) 


E^ - 


-\{-^,Dt2+2Ct2+DtJ, 


(59) 


F = - 


-\{Bt,+Dt2). 


(60) 



This exercise explicitly illustrates that two of the five in- 
troduced functions arc sufficient to fully characterize the 
polarization tensor. In particular, the latter is known en- 
tirely once the off-diagonal components (parametrized by 
77-^1 and 77^2) are known. 

In addition to the terms listed above, there are terms 
that are not constrained by the Ward identities, for in- 
stance 

7roo(0,q) = x(g), (61) 

which can be interpreted as a static susceptibility. When 
the time extent /3 is infinite, we have xio) ~ Q^ ^^ small 
momentum, due to the presence of a mass gap in QCD, 
and this guarantees that the photon remains massless. 
When /3 is finite however, x(0) is the thermal electric 
charge susceptibility, which does not vanish. Our parametriza- 
tion so far does not allow for this effect, and we must thus 
add the terms 



tM 



n^,' D xSg,o, IfiT D xS, 



9,0 



(62) 



where x ^i^d x ^-re 



to the generic expressions (40 44|p| 
two new coefficients. 

When some of the dimensions become infinite, non- 
analyticities in the momenta can appear. A familiar case 
is the limit where the spatial dimensions become infinite, 
but /3 remains finite. Then at small momenta, 77oo(9) ^ 

\\+d' 2 due to the diffusion pole (7? is the electric charge 
diffusion constant, see j30j for a review). 



At finite q, the coefficients B and B already capture the 
effect of the susceptibility x(q)- 
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In infinite volume, the vacuum polarization n{Q^) is 
uniquely defined, due to the constrained Lorentz structure 
of the polarization tensor. On the (5 x L^ torus however, 
we see that there are two obvious definitions of -/7(0), one 
that one may extract from the Ti representation, and one 
from the T2 representation. Indeed, to quadratic order the 
various components read 



7Toi(g) = Ar^qaqi, 
ni2{q) = AT^qiq2, 



(63) 

(64) 



^11(9) = (X + X)5qfl + Ar.ql - [Ar.ql + ^t,?'), (65) 



nooiq)^ ix-5x)Sq,o~MQ^ 



(66) 



Disregarding the susceptibility terms, they thus have the 
form dol) expected in infinite volume, except for TTn if 



qf ^ q^. In general however At^ ^ At^, and one has to 
specify how the extrapolation to g = is done. The vac- 
uum polarization extracted from the component i7oo((9) is 
the same as the one extracted from iloi, but one must be 
careful with Tin, given that the tensor structure appear- 



ing in Eq. ( 65 1 is not the one expected in infinite volume 
unless At-i = At2 ■ 

The equalities (451 to (60) can be used to test com- 



is thus convergenirl The Ward identity (67 1 can then be 
written 

/CXD 7 
— Ap33{io,k3,/3,L)=0. (69) 

-00 w 

For ^3 7^ 0, this sum rule constrains the finite-volume 
alterations of the spectral density relative to the infinite- 
volume situation. It may be useful in the finite-temperature 
context, where one wants to determine the distribution of 
the spectral weight Z\p33(w, fca, /3, L), see [30] and refer- 
ences therein. If one takes the spatial volume to infinity, 
^3 can be made as small as desired and the sum rule for 
Z\p33 then also holds for k^ — 0. We note that the sum rule 
(69 1 is simpler in QCD than in the Af = A super Yang- 
MiUs theory, because in the latter case the contribution of 
the scalar fields to the current contains a derivative, which 
means that the current itself is not invariant under a lo- 
cal symmetry transformation, and this leads to a contact 
term in the current correlator 1311. 



4 Computational strategies for Aa{Q'^) and 

^HLO 



puter programs. Once these tests are passed, a number of 
checks for finite- volume effects can be made. Apart from 

Ati = At2 , one can check for instance whether Bt^ van- 
ishes and C'ti = Dti ■, as one expects in infinite volume, 
and similarly in the T2 channel. 

Finally we note that if one calculates the isospin cur- 
rent correlator and introduces twisted boundary condi- 
tions for the quarks fields (see for instance [6 and Refs. 
therein), then the cubic symmetry is in general broken 
down even further, and the analysis above should be adapted 
accordingly. 



In lattice QCD it is customary to work with correlation 
functions which are functions of Euclidean time t and spa- 
tial momentum fe. This representation has the advantage 
that the low-lying states dominate exponentially at large 
t. Here we first consider the situation in infinite volume, 
deriving the relation between the Euclidean correlator and 
the spectral function p{q^) introduced in Eq. (pi). In the 
next section we show that a^ can be calculatea in terms 
of the (i, fe)-dependent correlator using a different kerner] 
We propose a way of treating different intervals of t differ- 
ently when computing a^ , combining the Euclidean-space 
calculation with the low-energy part of the spectral func- 
tion. 



3.2 A sum rule for the finite-volume spectral function 



The Ward identities k^II^y 



imply 



k3n33ik) = 0, fc = (0,0,0, fc3)- 



(67) 



This relation holds both in finite and in infinite volume. 
We now want to turn it into a sum rule for the spectral 
function P33 (introduced in Eq. I2| via the dispersion re- 
lation for the Euclidean correlator. The spectral function 
associated with the spatial current correlator 7T33 grows 
as fcp at large frequencies. Therefore a subtraction is nec- 
essary to obtain a convergent sum rule. Taking the dif- 
ference between the finite-volume and the infinite volume 
correlator leads to a subtracted spectral function oc kg^ 
at large frequencies, due to the absence of operators of 
dimension less than 4 in the operator product expansion 
of TT^t^. If /\/033 is the difference between the finite- volume 
spectral function and the infinite- volume one, the spectral 
representation 



4.1 The pion form factor in the timelike region 

The pion form factor in the timelike region can in principle 
be calculated based on Eq. (It]) in the threshold region, 



2m^± < \/s < 2 {ni^± + rrij^o ) . 



(70) 



Via Eq. ([5]), it can directly be compared to the experimen- 
tal results for the R ratio (for a precision comparison, care 
must be taken of QED corrections). The quality of the ex- 
perimental data can be viewed in Fig. (3) of Ref. [5T]. A 
further constraint on this region comes from the fact that 
the form factor F-,^{s) admits a convergent Taylor expan- 
sion in s with a radius given by the parameters of the p 
meson, (rn-p, Fp). The first two terms are well constrained 
by the pion form factor in the spacelike region. 



FAs) = 1 + ItI-s + 0{s^). 



(71) 



nsaih, /3, i) -^33(^3, 00, 00) 



dcu 



Ap33{uj,k3,/3,L) 



* We have written U and p as functions of the momentum 
variables that do not vanish. 

^ H.M. is indebted to Andreas Jiittner for bringing to his 
attention coordinate space methods for the calculation of a^. 
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This means that calculating the pion radius in lattice QCD 
also helps to constrain/check the spectral function in the 
threshold region. Although the region ( 70 ) only accounts 
for about 10% of the total a^^'^ , it is still important to 
control its contribution at the level of a few percent in 
preparation for the upcoming (g — 2)^ experiment at Fer- 
milab "52^. 

This program would be a relatively modest but cer- 
tainly valuable contribution of lattice QCD to the preci- 
sion determination of Aa{M'^) and a)^^'^ ■ 



4.2 Mixed-representation correlator 

We consider the positive-definite correlatoijj 

Git)^ f dx{jr{t,x)jf^\Q)). (72) 



correlator t^G{t) pictured in Fig. (16). Beyond say 2fm, it 
falls off rapidly to zero, asymptotically as e"^"*"*. We note 
that configuration-space Euclidean correlators have been 
used before to confront instanton [3 3) and holographic 
models |34j with experimental data via dispersion rela- 
tions. 

One can invert the Fourier transform (|73[), and ex- 



press the vacuum polarization ( 74 1 through the mixed- 
representation correlator, 



i7(c.2) = 



1 



Ate-"^'G{t). 



(79) 



At small cj, the vacuum polarization behaves as 

^ / dtG(t)-- / dteG{t) + ... (80) 



7T(c.2) ";^o 



We will now derive a spectral representation for G{t) in 
terms of the spectral function p(g^) of Eq. (3|. To this 
end, we note that the correlator G{t) can also be obtained 
from n^u{q) by Fourier transformation. 



The —2 term vanishes, as the corresponding integral rep- 
resents the quark number susceptibility of the vacuum (if 
one thinks of z as the 'time' direction). Thus 



7T(a;2) - 77(0) 



/OO J 
-^i7,,(c^,fc = 0), 

Now, the tensor structure ^ implies 



and secondly we can substitute the dispersion relation \10\ 
into Eq. (73 1. Noting that the 77(0) term only contributes 



dtG{t) 



(81) 



(73) 



(74) 



This integral is UV-finite by power counting (leaving in 
the term J dtG{t) = makes this explicit). We also note 
that G(i) is an even function of i, so that the expression 
is real, and we can write 



77(^2) -77(0) = 2 / AtG{t) 



1 — cos ujt 



(82) 



for t = 0, we obtain for i 7^ 

'•°° do; 



G{t) 



27r 



-UJ 



As 



Pis) 







s{s + LO^) 



(75) 



The integral is easily carried out and one obtains (again 
for t ^ 0), 



dtG{t) ujH^ - Asm^^ujt) .(83) 



The leading hadronic contribution to the anomalous 
magnetic moment of the muon can be ex p ress ed t hrough 
the vacuum polarization n{Q^), Eq. ([l2| and 14 1. Using 



relation ( 83 1 , we can write a}^^^ as 



1 r^ 
G{t) = -j As y5p(s)e-^l*l 



dcjoj p(w )e 



M 



(76) 
(77) 






Kit) 



mi_,t^ Jo 



00 

Att^ G{t) K{t) 

°° Au 



(84) 



OJ 



KEiio^)[LoV-Asm^f)].i85) 



Eq. (77) is the sought after spectral representation. One 
can of course derive the spectral representation of the cor- 
relator with non- vanishing spatial momentum, 



da;e-*-(j,(i,a;)jt(0)) (78) 

Auj (w^ - fc2 + fc2) p(u^ - k^)e-^W . 
|fc| 



Using the same parame triza tion of the phenomenological 
7?(s) ratio as in section (2.1), we obtain the dimensionless 



The kernel K{t) is dimensionless, proportional to t at 
small t and to 1/t at large t. The factor t^ was chosen 
because t^ G{t) is dimensionless and finite as t — >■ 0. The 
integrand of Eq. (84) is displayed in Fig. (|4]). 



4.3 Euclidean correlator vs. the R ratio 

We have already confronted the vacuum polarization cal- 
culated on the latti ce t o the R{s) ratio via the dispersion 
relation in section (2.1). Here we propose to do the same 



^ In Minkowski space, the field operators j^ are hermitian, 
jji ~ Jm- In Euclidean space, while Jq — jo, the spatial current 
is antihermitian jl = —jz- 



in the mixed (i, k) representation. This has the advan- 
tage that the correlator involves only on-shell states. In 
addition, the continuum limit is approached with O(a^) 
corrections (assuming the vector current is improved). 



David Bernecker, Harvey B. Meyer: Vector Correlators in Lattice QCD: methods and applications 




Fig. 3: The kernel of the spectral representation (77 1 



for three different Euclidean times, and the experimental 
i?(w2) ratio up to w = 2GeV Tf. 




in infinite volume, ob- 



Fig. 4: The integrand of Eq. 

tained from the phenomenological parametrization ( 18 ) of 

the R{s) ratio. 



with 



Consider the correlator G(t\ it is related to the spec- 
tral function p{uj'^) via Eq. (77 1. Depending on the value 
of t, the Euclidean correlator G{t) is sensitive to differ- 
ent energy intervals Aim. This simple fact is illustrated in 
Fig. (Is]). Confronting in this way the value of G{t) cal- 
culated on the lattice with the phenomenological value 
offers a way for the lattice practitioner to test the validity 
of the R{s) parametrization. Of course this test will only 
be useful if the correlator, as well as its statistical and sys- 
tematic error, are accurately determined. This also means 
that the chiral extrapolation, if any, must be under good 
control. Choosing t between 0.5 and 1.5fm yields sensitiv- 
ity to the region most relevant to the muon anomalous 
magnetic moment, encompassing in particular the p, w 
and (j) resonances. The comparison can therefore provide 
a useful test of the claimed accuracy of the phenomeno- 
logical approach. 



a<(to,L) = Aa^m^ / dtt^ G{t) K{t) (87) 



Aa^{to,ti,L) = Aa^m^ I dtr G{t) K{t). (88) 

Jto 
/•OO 

a>{ti,L)=4a^m^ dtt^ G{t) K{t). (89) 

Jti 



The short-distance contribution a1^ (to) — lim^ 



'iL,t) 



,HLO 



can be calculated in perturbation theory; the perturbative 
series is known to have good convergence properties in the 
vector channel. For the purpose of calculating a^, the prac- 
tical question will be whether there is a choice of to where 
the perturbative series shows good convergence and the 
discretization errors are smalQ 

A state of energy E makes a contribution of order 
e~'^*i to ar^(ii, L). With a solid understanding of the low- 
energy spectrum, one can therefore analyze this contribu- 
tion for large enough ii and in particular its finite-volume 
effects. The larger ti, the stronger low-energy states dom- 
inate, but on the other hand the numerical importance of 
this contribution is reduced. Therefore one would like to 
choose ti as small as possible while the contribution is still 
dominated by the analytically tractable low-lying states. 

In appendix ( [A| we show that the finite-size correction 
on the contribution coming from p{s < 4(171,^^ +171^^^)'^) to 
the long-distance part of the correlator is large (by long- 

distance, we mean t ^ (2m7r)^^ and t ^ L/tt). Figure 

The most relevant ones are the lepton mass m^, the ^ illustrates that this contribution represents 50% of the 

full correlator starting at i « 2.9fm. The contribution of 
G{t > 2.9fm) to a^^° is modest (see Fig.li]), but certainly 
not negligible if one aims at a precision at the percent level 



4.4 Combining spacelike and timelike correlators 



There are many scales entering the integral ( 84 ) yielding 



pion mass TOtt, the box size L and the integration variable 
t. An advantage of working in the mixed (t, k) representa- 
tion is that a Hamiltonian interpretation of the expression 
is straightforward. This motivates us to separate the con- 
tribution of different t- intervals to the integral ( 84 1 in the 



following way. We write the integral as a sum of three 
terms, 



or better on a^^°. 

In view of these finite-size corrections, it appears prefer- 
able to treat the contribution a^(ii) differently, if ii is 



HLO 



^ An analogous issue arises when computing a^ in four- 
momentum space, where perturbation theory is found to match 



a1^(tQ, L) + Aafj_{tQ,ti, L) + al^[ti, L), (86) the lattice vacuum polarization down to about 3GeV^ [4]. 
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Fig. 5: The fraction of the q — Euclidean current 
correlator G(t) coming from the region 277i^± < -^s < 
2(77i7r± + ni^o) of the spectral function (in infinite spatial 
volume). See Eq. (JTh. 



chosen large. Eq. ( 121 1 is an exact spectral representation 



of this contribution. The infinite- volume spectral function 
p{oj^) can be calculated on the lattice below the four-pion 
threshold via Eq. ([6][7|. Therefore, if the states below the 
four-pion thresholdsaturate the current correlator beyond 
ii, a^l^i) can be obtained from first principles through 



Eq. (1211, while Aa^ is calculated by integrating the Eu- 
clidean correlator as in (88). This is an example of the 



a different diniensionful quantity to 'set the scale'. It is 
however not easy to come up with a quantity which is 
both accurately calculable on the lattice and accurately 
extracted from experimenl|^ A quantity that has proved 
very useful is the Sommer reference scale ro, defined from 
the static quark potential [3S]. While it is accurately cal- 
culable on the lattice, its value in the real world is not 
known precisely, since its extraction from the upsilon spec- 
trum requires introducing a potential model. This proce- 
dure leads to a certain degree of ambiguity (up to ^ 8%, 
to be conservative). 

Here we propose a new reference scale tq based on the 
vector current correlator, which we believe satisfies the re- 
quirement of being accurately calculable. Given the level 
of experimental effort that has gone into the measurement 
of the R{s) ratio particularly in the past decade, we be- 
lieve that the value of tq can also accurately be extracted 
from the compiled data of e^e~ annihilation experiments. 
Figure ([6]) displays the electromagnetic current correlator 
G{t) obtained from our simple parametrization of the R 
ratio, Eq. (pi. 



potentially powerful interplay of a correlator in the space- 
like and in the timelike region. The difhculty with this 
strategy is that numerically, the saturation at (say) the 
95% level by the states below the four-pion threshold only 
begins around 5.3fm, at which point the contribution to 
aj)""° is already extremely small. Choosing ii as large as 
5.3fm is therefore not a very attractive option. The strat- 
egy of splitting up a^ into three contributions and us- 
ing the spectral representation for the long-distance part 
would become more attractive if one could determine the 
spectral function up to somewhat higher energies, but the 
formalism has not yet been developed to handle the mix- 
ing of two-pion with four-pion states. 

The analysis outlined above shows that the relative 
finite-size effects on the long-distance contributions to the 
vacuum polarization and to a^^'^ are large. Fortunately, 
the size of these contributions is small (albeit non-negligible) 
compared to the total a}^^^ . The finite-size effects asso- 
ciated with the higher-lying states remains however un- 
known, and it should be investigated. 



5 A new reference scale for lattice QCD 

When making predictions for hadronic observables in lat- 
tice QCD, a mandatory step is to calibrate the length of 
the lattice spacing a in physical units (fm). In principle, 
calculating the proton mass Mp in lattice units, and equat- 
ing it to 938.272MeV yields the desired value. However, 
the proton mass is difficult to calculate accurately on the 
lattice at light quark masses, and in practice one chooses 



At long time separations (somewhere beyond Ifm), the 
vector channel correlator becomes noisy in Monte-Carlo 
simulations. At short distances, one is confronted with 
cutoff effects from the lattice, and secondly, QCD is ap- 
proximately scale-invariant, so that the sensitivity to the 
confinement scale is low. This dictates that one should 
choose a reference time-scale somewhere between 0.5fm 
and Ifm. 

Another computational aspect is that it is significantly 
easier to calculate correlators in channels where no Wick- 
disconnected diagrams appear. With the commonly used 
statistical sampling methods, connected diagrams exhibit 
better signal-to-noise ratios. The correlator of the electro- 
magnetic current will contain the difference of the strange- 
quark disconnected loops and the light-quark disconnected 
loop. Therefore, for the purpose of statistical accuracy it 
is preferable to use the isospin current, in which discon- 
nected diagrams cancel out (assuming exact isospin sym- 
metry). On the phenomenological side, this means that 
one should only include those final states with isospin 
/ = 1 in the evaluation of the i?-ratio. We therefore define 



Ri[s) 



a{e'^e 



hadrons I / = 1) 



47ra(s)2/(3s) 



(90) 



In particular, the bulk of the uj contribution, which decays 
into an / = three-pion state, should not be included; nor 
should the (f) meson be included. The bulk of the low- 
energy isospin-current spectral function will thus come 
from the two-pion channel. The restriction to / = 1 final 
states can be done in a model-independent way as long 
as isospin breaking effects in QCD are negligible to the 
desired degree of accuracy. The number of open channels 
between 1.4GeV and 2.0GeV increases rapidly, and above 
2GeV, practically only inclusive measurements are made. 
Therefore it is preferable to suppress the high-energy re- 
gion by choosing a relatively large tq ■ 



If the goal is only to calibrate the relative size of two lattice 
spacings, the second requirement is not mandatory. 
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Fig. 6: Top: the correlator of the electromagnetic current 
J™ at zero spatial momentum in infinite volume obtained 
from the phenomenological paramctrization of the R{s) 
ratio. Bottom: the same for its isospin component jj,. 



We thus define the Euclidean correlator Gi{t) analo- 
gously to G{t) (Eq. 72 ), but replacing the electromagnetic 
current j°™ by its isospin component 



ji = \{u'y^,u-d-if,d). 
We define the time scale tq by the equation 



dcf 



To^Gi(to) = 0.027 



(91) 



(92) 



because it leads to tq in the desired range 0.5fm to l.Ofm. 
Figure (l6| displays the electromagnetic current corre- 
lator G{t) obtained from our simple paramctrization of 
the R ratio, Eq. (18 1, as well as the corresponding isospin 
current correlator Gi{t) based on the model 



i?i(s) = e{^/s - 2m^±) 6l(4.4m^± - ^/s) 



(93) 



1- 



4m^ 



3/2 



(0.6473 + /o(ys)) 



+0(^5 - 4.4m,±)0(M3 - ys)/i(\/^) 
+3((i)2 + (i)2)0(V.-Af3). 

(the values used for the parameters are those in table fTl). 
Based on this approximate form of the spectral function, 
we find 

To ~ 0.85fm. (94) 



We postpone a more accurate phenomenological evalua- 
tion of To and its uncertainty to a future publication. At 
a practical level, a nice feature of the tq definition is that 
no derivative must be taken and no fit need be performed, 
one will 'merely' have to perform an interpolation to the 
point Tq. 

Of course, many reference scales can be defined in a 
similar way. As a slight variation of the definition ( 92 ), one 



could also use the time component of the current with a 
non- vanishing spatial momentum fe. The advantage is that 
the correlator does not vary as fast at short distances, only 
as k It instead of \/i^ . The spectral representation then 
readsQ 



da;e-^'=-^(jo(i,a;)jo(0)) (95) 

/•DO 

J\k\ 



A reference scale based on this correlator might be an at- 
tractive alternative depending on its signal-to-noise ratio. 
If one manages to calculate the disconnected diagrams 
accurately, then a definition based on the electromagnetic 
current could become more attractive. For instanc e, de fin- 



dcf 



ing r^j^G(rcm) — 0.0336 leads again, using Eq. (18), to 



Tern ~ 0.85. This has the advantage that the selection of 
unit-isospin final states in the R[s) ratio is not required, 
nor is the assumption of exact isospin symmetry. 

Finally, we remark that as with any another reference 
scale, a chiral extrapolation is necessary unless a simula- 
tion is performed at physical quark masses. We leave this 
question for future investigation, but note that the prop- 
erties of the p meson probably play an important role in 
this respect. 



6 LiJscher formula with twisted boundary 
conditions 

For the evaluation of formula (ItI) the scattering phase shift 
(5i of the 1 = 1 = 1 channel must be known as a func- 
tion to obtain its derivative. Liischer showed precisely |16| 
how individual values of the scattering phase can be re- 
constructed out of the two particle energy spectrum in- 
side a finite volume. The problem that arises is that in 
practice, for a given volume size, only a few values of the 
scattering phase shift can be reconstructed. If the scatter- 
ing phase shift is needed as a function of momentum the 
computational cost increases rapidly since the results of 
many simulations with different volume have to be com- 
bined. Rummukainen and Gottlieb [36] showed that by 
extending the formalism to non vanishing center of mass 
momenta more points of the scattering phase shift can 
be extracted per volume. A further generalization is the 
use of twisted boundary conditions that introduce a new 
continuous parameter, the twist angle, to the simulations 



It is easily obtained from Eq. ( 78 1 by using the Ward iden- 
tity ^^J,jf_, = 0. 
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(see for instance [37]). The twist angle can then be used 
to modify the momentum almost continuously and there- 
fore enables one to "scan" the scattering phase shifts. An 
elegant proof of the extended version of the Liischer for- 
mula was already given in [T^ . Our derivation closely fol- 
lows the original one in [TS] and is therefore only sketched 
here. First we introduce the wave function in infinite vol- 
ume. It describes two bosons of equal mass, and obeys the 
Schrodinger equation 



r»,2 



Pi 

2fj, 



V{\r\) 



!f(r) = E^{r) 



(96) 



where the potential V^(|t|) is spherically symmetric and 
has a finite range R. The wave function can be extended 
into spherical harmonics with its radial components satis- 
fying the radial Schrodinger equation 



dr'^ 



2 d 

r dr 



1(1 + 1) 



k^ - 2^V{r 



^imir) = 0. (97) 



In the region where the potential vanishes, the solution 
is given by ^im{r) = ai{k)ji{kr) + j3i(k)ni[kr), where the 
constants a; and /?/ determine the scattering phase shifts 
via 

2ts,{k) _ aijk) + if3i{k) 



ai{k) - iPiik) 



(98) 



We now enclose both particles in a box of size L x L x L 
and impose twisted boundary conditions 



W{r + nL) =e*"'^if(r). 



(99) 



where the triplet cf) of twist angles was introduced. The 
potential must also be replaced by a periodic version 



ydr)^ J2 Vi\r + nL\). 



(100) 



ngZ 



Inside the box the energy spectrum is now discrete, with 
the energy- momentum relation still given by -B = fc^/2/i. 
As in the periodic case, the Schrodinger equation in the 
region of vanishing potential f2 reduces to the Helmholtz 
equation, 

{A + k^)^{r)^0. (101) 

In this outer region i? the eigenfunctions of the Hamilto- 
nian now have to be solutions of the Helmholtz equation, 
that are expandable in spherical harmonics and have ra- 
dial components 



I'lm (r) = bira («( {k)ji (kr) + Pi {k)ni [kr)). 



(102) 



According to a theorem by Liischer (proof given in ap- 
pendix A of [16j ) for each solution of the Helmholtz equa- 
tion in fi that can be expanded in this way there exists 
a unique eigenfunction of the Hamiltonian that matches 
this solution in fi. Finding the general solution of the 
Helmholtz equation therefore suffices for determining the 
formula for the scattering phase shifts. 

To simplify the derivation a bit, the momenta of the so- 
lutions of the Helmholtz equation are now assumed not to 



belong to the singular set Fg — {k £'&\k = ±'^\n\, ior n <E 
1? } . Momenta in this set would allow plane waves as solu- 
tions of the Helmholtz equation, which would complicate 
the derivation a bit. For the following steps an angular 
momentum cutoff A is introduced, so that only the par- 
tial waves with an angular momentum smaller or equal to 
A feel the presence of the potential. 

The solutions we are about to construct should sat- 
isfy two conditions. First they have to satisfy the twisted 
boundary conditions ( 99 1, and second they should be bounded 
by a power of r~^ near the origin 



lim |r^+^tZ/(r)| < oo. 

r— fO 



(103) 



The ansatz for finding the general solution is now the 
Greens function 



G^(r;fc2)=L- 



per* 



where the momenta are elements of F'^' = {p G 



(104) 



P = 



It: 



n + L^^<p, n e Z'^}. That this function satisfies the 
twisted boundary conditions can be shown straightfor- 
wardly, and it is also a solution of the Helmholtz equa- 
tion for r 7^ (mod L) . By comparison with the spherical 
Bessel functions one finds the behaviour of G'^{r; fc^) near 
the origin to be 



G"^(r;fc2) = — no(fcr) + G'^(r;A:2) (105) 

47r 

where G'^{r; fc^) is the regular part. Further solutions may 
be constructed by using the function yim{x) — r^Yim{d, (t>) 
to obtain the derivatives 



GUr;k^)^yU'^)GHr;k'). 



(106) 



One can then show that the Gf^{r; fc^) form a complete, 
normal basis of the solutions of the Helmholtz equation. 
Therefore the solutions that were searched for can be con- 
structed as a linear combination of the Gf^{r; fc^) 



A I 



•^W=E E vimGUr;^). 



(107) 



1=0 m=-l 



Like in Liischer's paper these solutions now have to be 
expanded in spherical harmonics and this form has to be 
compared to the expansion containing the spherical Bessel 
functions. The expansion is similar to the one with peri- 
odic boundary conditions 



GLir;k') 



(-l)'A:'+i 



{niikr)YiUO,cl>) 



An 

oo I' 

E E <nMm'{'lik))jl'{kr}Y,^'{eA)}, 

(108) 



l'=Om'=l' 
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where q 



^ and Mt,,. 



I ('?(^)) is given by 



M 



Irn.l'rrf 



(9) 



4)' 



r3/2 



E 



J 



A^ ffi+i 



Ci 



]=\i-i'\ ^=-'J 



qj 



lm,js,l'm 



'2t(l;g^). 

(109) 

While the coefficients C'im,js,i'm' are given in fl6J and can 
easily be calculated for a given set of indices, the Zeta- 
functions Zf^^{l;q'^) (defined below) need to be deter- 
mined numerically. From the comparison of the expan- 
sions one finds (for the choice cp — {(f), (p, (p) so as not to 
break the cubic symmetry) the usual Liischer formula 



det 



^2iS 



M'i'in) 



M'^(7^) - 







(110) 



where now the M*^ are determined by a different Zeta- 
function. This formula is already the reduced one that was 
obtained by projecting on one of the irreducible represen- 
tations {TZ ^ Af, A^, E^, T^, T^) of the cubic group, 
as is described in [161 . The remaining task now is to cal- 
culate the Zeta-functions Zf„_,{s]q^) that are defined by 

(111) 



""hn 



for Res > 1, with /"^ = {r e M?'\r = n+ (27r)-V, " e 
Z^}. Like in Liischers original work the Zeta-functions are 
best calculated numerically using integral representations 
of the form 

|r|<A 
re/"*" 



+ (2-)^ f 
Jo 



dt 



e''nt\t,o) 



<t>,\. 



5m6.. 



W-'mO 



(47r)2i3/2 

(112) 



Here the reduced heat-kernel nl^ (t, x) of the Laplace 
operator on a torus with twisted boundary conditions of 
size L — 2n appears, and A must be chosen so that A^ > 
Re q^ . Depending on the value of t the heat kernel has two 
different representations based on 



H^it,x)^{4Trtf/^ 



n& 



e^'^-'t'e 



i(x-27rn) 



(2^) 



-3 \ ^ ir-x — tr 



(113) 



the first converges for t < 1 and the second one for t > 1. 
With those two representations the reduced version of the 
heat-kernel that is needed for the integral representation 
of the Zeta-function can then be defined by 



,4,,\ 



nt^{t,x) = {-iyyUV) nnt,x) 



(2^) 



— 3 \ ^ ir-x—tr 
\r\<X 




(a) periodic boundary conditions 




(b) twisted boundary conditions (L = 20) 

Fig. 7: The reconstructed scattering phase shift 6i{k) us- 
ing the Liischer formula with twisted boundary conditions. 
For futher information on this example see appendix iBl 



Using for instance the program Mathematica the Zeta- 
function can then be calculated numerically. 

Figure [7] illustrates the advantage of using the twisted 
boundary conditions. For a simple quantum mechanical 
example with a potential well (see appendix mn the re- 
constructed scattering phase shift Si{k) is shown. Fig. 7a 
shows the scattering phase shift reconstructed using peri- 
odic boundary conditions and several volumes, while for 



fig. 7b twisted boundary conditions were used. In the sec- 
ond figure the volume was set to be L = 20, while the twist 
angle ranges between </> = and = 0.4. In both figures 



the blue curve is the scattering phase shift calculated in 
infinite volume that was added as a comparison to the 
numerically calculated data. This simple example already 
shows the advantages of the twisted boundary conditions, 
since the number of points that could be reconstructed 
from one single volume could be improved. Also the dis- 
tances between the individual points is much smaller, so 
that is would also be possible to calculate the derivative 
of the scattering phase shift, which is needed for equation 
(114) 0. 
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7 Conclusion 

Given the highly accurate determinations of the R ratio in 
e+e~ experiments, where now even the dominant isospin 
breaking effects are taken into account, it is a challenge for 
lattice QCD to have a phenomenological impact on the de- 
termination of the hadronic vacuum polarization. In view 
of the upcoming (g — 2)^ experiment at Fermilab |32j . an 
accuracy of well below 1% is called for on a^^'~'. We have 
outlined two alternatives to the direct calculation of the 
vacuum polarization. One option is the direct determina- 
tion of the spectral function p{s) in the threshold region 
along the lines of [T^ . We described the determination of 
the scattering phases in the vector channel using twisted 
boundary conditions in section (rol. A second promising 
approach is the calculation of the vector correlator G{t) 
(defined in Eq. 72) in the region 0.5fm to 1.5fm or so. 



which provides a check on the R ratio in the energy re- 
gion that makes the largest contribution to a^^'~' . We also 
considered the option of integrating G{t) from Iq to ti, 
with the short-distance contribution treated in perturba- 
tion theory and the long-distance contribution obtained 
from the pion timelike form factor via the spectral repre- 
sentation, but it appears that the spectral function would 
have to be determined up to higher energies than the four- 
pion threshold for this strategy to be practical. 

As a new idea we proposed to turn the table around 
and exploit the accurate knowledge of the vector spectral 
function to define a reference scale tq in QCD which can 
be determined accurately and reliably. We believe that 
the only real difficulty could arise from finite-size effects, 
which should therefore be investigated numerically and 
analytically. 

We have made first steps to study the finite-size ef- 
fects affecting the calculation of the vacuum polarization 
and the hadronic contribution to {g — 2)^. In section (pi), 
we have analyzed the tensor structure of the polarization 
tensor on the torus, and find that for generic values of the 
momentum, two functions invariant under the transforma- 
tions of the cubic group H{3) characterize the vector cur- 
rent correlator. This allows one to perform various tests for 
finite-size effects. Based on current conservation, we also 
derived a sum rule for the difference of the infinite- volume 
spectral function and the finite-temperature/finite- volume 
spectral function, Eq. (69). Then focussing on the spatial 



current correlator projected onto zero spatial momentum, 
we analyzed the contributions of different states as a func- 
tion of Euclidean time. At large times, the relative finite 
volume effects become order unity, due to the discrete- 
ness of the two-pion states on the torus. This contribution 
however only represents a few percent of a^ , due to the 
dominance of the p resonance in this quantity. We plan to 
study the finite-size effects of this contribution in the near 
future. 
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A Finite-volume effects on the low-energy 
contribution to a^^^ 

In dealing with the different scales of the problem we will 

consider that 

•III — 

= 0(1). (115) 



rrhn 
m,, 



Because of this relation, the kernel K{t) cannot be ex- 
panded, neither in a small-i expansion, nor in a large-t 
expansion. 

In the chiral regime, the lowest-lying states contribut- 
ing are two-pion states. Let us call the ground state energy 
Eq . Roughly at Eq + 2m.^ begin the states with a four-pion 
component in their wavefunction which is not exponen- 
tially small in the volume. Since it is not known how to 
relate the relative weight of these components to inelastic 
TTTT scattering in infinite volume, we choose ti in such a 
way that the contribution of these states is suppressed, 

2m^ti > 1. (116) 

We assume that the box is large, in the sense that 



If we also choose ii such that 



ii > 



(117) 



(118) 



then only the ground state contributes to G{t) for t > 
ti- Since all factors in the integrand are slowly varying, 
except for the exponential e~^*, the contribution to a^ is 
approximately 



a>iL,ti) ~ Wm^tlK{ti)\Ao{L)\' 



-Ecjti 



En 



(119) 



where Aq — L^/^(vac|j(a;)|-0o) is the matrix element of the 
vector current between the lowest-lying (two-pion) state 
in the box and the vacuum; it is related to the pion form 
factor through Eq. (pi. In infinite volume, the situation is 
different, because there are states arbitrarily close to 2mT^, 
and we must integrate over them. We define 



au(^i) = lim a>{L,ti). 



(120) 



Using the spectral representation of the correlator, one 
finds 

/■OO /"OO 

a>{ti)^Aa^m^l dujuj^p{uj^) dtt^ K{t) e''^* (121) 



4:amf_,K{ti)fl dcj ujp{uj)e 
Jo 



2\ —ujti 



(122) 
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The large value of ii dictates that only the small w region 
contributes, where the spectral density coincides with the 
timelike pion factor, see Eq. (|6|. One then finds 



to calculate the spectrum in the box. The general solu- 
tion of the Helmholtz equation is used to define a solution 
piecewise inside the box. 



:{ti)^Aa^m^,tlK{ti) 



|F,(2m,)|2 3^ 



-2mTrti 



487r2 



1/2,5/2 ■ 



^ ^ ^ ^ (123) 

Comparing Eq. (119 ) and ( 123 ), we see that the expression 



obtained in finite volume is parametrically different, be- 
cause exp(— ii^o(-^)^i) ^ exp(— 2TO7rii). The finite-volume 
effect is order unity for this long-distance contribution. 
If the box is very large, rriTrL ^ tt, the contribution 



Z\a^(L,ti) for ti 



< L 



(but still 2m^ii > 1) will be af- 



V{\r\) = -a-0{R^\r\), 



(124) 



where R is the range of the potential and a its strength. 
For a comparison to the scattering phase shifts calculated 
using the Liischer formula the scattering phase shift in 
infinite volume is needed. To solution to the Schrodinger 
equation 



1 
2/i 



A + V{\r\) 



•^(^) - ^"^C^) 



(125) 



is best defined piecewise 



if(r) 



CX3 / 

EE 

/=0 m=-l 



Ylr 



\ aiji{kr) + I3ini{kr) r>R 
'\3i{kr) r<R. 

... ^^2^) 
Requiring smoothness of the solution and its derivative 

at r = i? leads to a system of equations for the coeffi- 
ciants a/(fc), /3i{k). Once these are calculated the scatter- 
ing phase shift can be calculated using the formula 



Si{k) — arctan 



(127) 






\r\>R 



\r\<R 



-E E -in.Gt{r-k% 

1=0 m=-l 
A I 

^ ^ ai7nYim{d,<P)ji{kr) 

1=0 m=-l 

oo / 

+ ^ ^ aimYim{0,(l))ji{kr). 

l=A+l m=-l 



(128) 



(129) 



fected only by a small relative finite- volume effect. In this 
regime, the energy gaps between the two-pion states is 
much smaller than myr, so that the sum over states is a 
good approximation to the integral over momenta in infi- 
nite volume. This regime was investigated in |38j . How- 
ever, this regime is hardly realistically achievable with 
present computing resources, because the condition ttIt^L ^ 
TT implies that L must be lOfm at the very least when m-^ 
is set to its physical value. 



B Illustration with a potential well 

In section ^ the scattering phase shifts calculated for 
the quantum mechanical potential well were shown as an 
example. This appendix describes how these scattering 
phase shifts were calculated. 

The potential well is of the form 



The advantage of this ansatz is that the general solution 
of the Helmholtz equation that is used already satisfies 
the twisted boundary conditions 



i^(r + ni) = e*"-'^tZ/(r), n e Z^ 



(130) 



To simplify the calculation and to study the effects of 
higher scattering phase shifts the angular momentum cut- 
off A was introduced in the same way as in section ([6]). 
Before matching both pieces of the solution at the bound- 
ary r = R, the solution outside the range of the potential 
can be rewritten as 



A I 



oo / 

E E 5^'™'W){ 



(131) 



l'=Om' = -l' 

1=0 m=~l 

+ e{A-i')vi. 



Vlr. 



(~)'fc'+l 

47r 



■M 



lm:l'' 



(132) 



(-)''A:''+i 



47r 



-ni 



(fcr)}. 



Since the Yim{0,(t>) are normal to each other one gets a 
homogeneous system of equations when equating the two 
pieces of the solution at r = i?. The important point to 
note is that the equations for Z > yl do not influence the 
spectrum, they merely determine the aiyA,m in terms of 
the vi<:A,m- The latter coefficients are determined by the 
equations I < A. One has yl -I- 1 homogeneous equations 
for as many variables. According to a well known theorem 
a non-trivial solution to a set of homogeneous equations 
only exists when the determinant of the corresponding ma- 
trix vanishes. This condition yields the allowed fc-values. 
For yl = 2 for instance, the determinant of the matrix of 
coefficients is of the block-diagonal form 



det[AO] • det[Al] • det[A2] = (133) 



AO 











Al 











A2 



In order to calculate the scattering phase shifts using so that it seperates. Since we are interested in the scatter- 
the formula derived in section ([6]) the first step now is ing phase shift Si we only need to consider the determinant 
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det[Al]. The system of equations that needs to be solved 
is 



= vi-i — [ji{kR)Mi-i.,i-i 
47r 



ni{kR)] +ai-iji{kR) 
(134) 



= Vio^[ji{kR)Mi„,io + ni{kR)] + a^a3i{kR) (135) 

^ = vii — [h{kR)Mii,ii+n^{kR)\ +anJi{kR}. (136) 

In this example the roots of the determinant were found 
by plotting the determinant as a function of k using Math- 
ematica and then using the function FindRoot. The calcu- 
lated k values can then be inserted into the Liischer for- 
mula for twisted boundary conditions derived in section 
(|6| to calculate the scattering phase shifts. The parame- 
ters that were used to obtain the figures [7a||7b| were a = 2, 
R = 2 and /i = 1. 
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